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ABSTRACT 

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio 
galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry 
from Spitzer, Herschel and LABOCA we fit the rest-frame 5 — 300 /im emission using a two component, 
starburst and active galactic nucleus (AGN), model. The total infrared (8 — 1000 /im) luminosity of 
this galaxy is (1. 97 ± 0.28) x 10 13 Lq with (1.17±0.27) and (0.79 ± 0.09) xl0 13 L Q due to the AGN and 
starburst components respectively. The high derived AGN accretion rate of ~ 20 % Eddington, and 
the measured star formation rate (SFR) of 1390 ± 150M Q yr _1 , suggest that this massive system is in 
a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger 
in a dense environment. The accretion rate is sufficient to power both the jets and the previously 
observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust 
its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with 
nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds 
of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, 
and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a 
radio-loud AGN phase. We also note tentative evidence for [CII]158/im emission. This paper marks 

the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of 

the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies. 

Subject headings: galaxies: active, formation, high redshift, individual (MRC 1138-262) 
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1. INTRODUCTION 

There is now a growing consensus that high redshift 
radio galaxies (HzRG; z > 1, £ 50 omhz > 10 26 WHz _1 ) 
are markers of the early formation of the most massive 
galaxies, located wit hin the highest peaks o f dark 
matter over-densities (jMilev fc De Breuck! |200§ ). Not 
only do HzRGs tend to reside in significantly over- 
dense environment s (e.g ., Ilvison et al.l 120001: iStern et al.l 
italT 1200' 



20031 iVenemans et al.1 120071: IGalametz et al.1 I201C 
Tanak a et al.l 120111: iGalametz et al.l 120121: iMavo et al.l 
2012D . but the mass of their stellar component 
also tends to be high as indicated by the remark- 
ably tight scatter between t he observed K— band 
magnitude and redshift (e.g.. iLillv fc Longairl Jl984t 
Eales et al.lfl997l Parvis et al.l 120011: iWillott et al.1 120031 : 
Rocca-Volmerange et al.l I2004t IBrvant et al.l 120091 ) . 
The high stellar masses have been confirmed more 
directly with rest-f r ame near-infrared ob servations 
(jSevmour et al.l 120071: iDe Breuck et al.l 12010ft . Recent 
observations of the high mid-infrared luminosities of 
HzRGs indicate that the black hole s in these sources 
likely have h i gh a c cretion rates (IQgle et al.l 120061 : 
iSevmour etldl 120071: iFernandes et al.l 120 111) . HzRGs 
are often found with evidence of high star formation 
rates (SFRs) from de t ections in the sub- millimetre (e.g., 
Archibald et al.ll200il; IStevens et al.ll2003t iReuland et al.1 



20041 ). from their extended Lyman a emis sion indicating 
large reservoirs of ionized ga s (e.g., IReuland et al.l 
120031: iVillar-Martm et al.l 12003). from their CO emis- 
sion indicating reservo i rs of molecular H 2 (e.g., 
iPapadopoulos et al.ll2000l: lEmonts et al.ll201lD. from hot 



star ultraviolet absorption lines (e.g., IDev et al.l fl997l) 
and f r om their mid-infr ared spectra (e.g.. ISevmour et al.1 
120081 : IQgle et aLll2012t ). Kinematic measurements from 
ionized gas have also shown that HzRGs have extremely 
massive central black holes (3 x 10 9 — 2 x 10 10 M Q - 
iNesvadba et al.ll2011l ) consistent with the observed high 
stellar masses and accretion rates. All of these obser- 
vations confirm the picture that HzRGs arc markers of 
massive galaxies co-forming with a central black hole, 
residing in peaks of dark matter over-density. 

The exact formation mechanism of massive galaxies, 
and how the frequently observed radio-loud phase is re- 
lated to black hole and galaxy growth, is often debated. 
A connection between the growth of the central black 
hole and the host galaxy has been proposed to recon- 
cile semi-analytical models with observations. A form 
of 'feedback', associated with a radio-loud phase, has 
been proposed to truncate the growth of the bulge and 
central black hole growth, thereby preventing the over- 
production of massive galaxi e s being formed in these sim- 
ulations (jBo^EiEaDlinQl ErotonetaDEnQl). This 
idea has gained popularity with the kinematic observa- 
tions suggesting powerful outflows from many HzRGs 
(e.g., [Nesvadba et al. 2006), although most models favor 
scenarios where the winds are radiatively driven. 

In order to test how phases of powerful radio-loud emis- 
sion influence, or are associated with, the formation of 
massive galaxies at high redshift, accurate measurements 
of the black hole accretion rate and the SFR are required. 
Due to varying, and difficult to determine, levels of ob- 
scuration at UV and optical wavelengths, as well as con- 
tamination from the active galactic nucleus (AGN), the 



best wavelength regime in which to obtain such mea- 
surements is the infrared (IR) where the spectral en- 
ergy distributions (SEDs) of such sources also tend to 
peak . With the Spitzer Space Telescope (jWerner et al.1 
2004) we conducted a systematic study with eight pho- 
tometric bands (3.6 — 160 /jm) of 71 HzRGs with red- 
shifts, z, in the range 1 < z < 5.2 in order to measure 
their near- to mid-IR rest-frame SEDs. This photom- 
etry allowed us to estimate separately the stellar and 
AGN contribution at these wavelengths, thereby con- 
firming th eir high stellar masses and AGN mid-IR lu - 
minosities (jSevmour et al.ll2007l : IDe Breuck et al.ll2010D . 
In order to measure the complete IR SED, we are now 
extending this study to far-IR wavelengths as the Her- 
schel Radio Galaxy Evolution Project (Projet He RGE) . 
The Herschel Space Telescop^ (jPilbratt et al.l I2010D 
has already demonstrated that far-IR measurements can 
relia bly be used to measure SFRs in powerful AGN 
(e.g^jH ardcas tle et alJl2010HHatziminaoglou et al.ll2010l : 



iShao et al.l 120101: ISevmour et al.l 1201 ID . By combining 
these data with those from Spitzer we can, for the first 
time, accurately measure the total bolometric luminosi- 
ties and far-IR contribution, and then estimate the rela- 
tive AGN and starburst IR luminosities in HzRGs. 

In this letter we present the first results from Pro- 
jet HeRGE on MRC 1138-262 at z = 2.156, otherwise 
known as the Spiderweb Galaxy due to the preponder- 
ance of satellite galaxies, 'flies', reve aled by deep Hub - 
ble Space Telescope (HST) imaging (Mi ley et alJl2006[ ). 
This HzRG has a powerful AGN as implied by its high 
X-ray lum inosity (i^-iokeV = 4 x 10 45 ergs _1 ~ 1.0 X 
10 12 Lq - iCarilli et al.l I2002T ) . its extr eme radio emis- 
sion (L 50 omhz = 1-2 x lO^WHz" 1 - IDe Breuck et all 
IMoft. and evidence of a high SFR ( ~ 1200 A^yr" 1 
- IStevens et~aTl 120031: IQgle et al.l |2012T). There i s also 
evidence of a powerful outflow (jNesvadba et al.l 120061 ) 
potentially driven by the massive central black hole 
(jNesvadba et al.ll2011[ ). A study of the proto-cluster en- 
vironment of the Spiderweb Galaxy shows that many 
members will be stripped by, or me rge with, the cen- 
tral galaxy by the present epoch (lHatch et al.l 120091: 
Dohc rtv et all MM iTanaka et all 120101 ). iKuiper et al l 
(j2011| ) suggest that the measured dynamical structure 
around the Spiderweb Galaxy can be interpreted as the 
merger of two massive haloes and that the combined sys- 
tem has not yet become viralised. In this paper we aim 
to measure more accurately the AGN accretion rate and 
SFR, the bolometrically dominant components, thereby 
obtaining an improved understanding of the evolutionary 
state of this HzRG. 

In Sj^jwe present the IR data used in this work and in <J3] 
we present the SED fitting of these data with an AGN 
and starburst component. We discuss our results, and 
the overall evolutionary state of the Spiderweb Galaxy 
in ||4] and conclude the paper in fj5j Throughout this 
paper we use the 'concordance' ACDM cosmology with 
H a = 70Mpc" 1 kms"\ Q A = 1 - Q M = 0.7. 

2. DATA 



22 Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 
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The data used in this work comprises observed in- 
frared photometry across 20 — 1000 /mi. At the shortest 
wavelength, 24 /am, t he photometry comes fro m Spitzer 
and was published in iDe Breuck et al.1 ([20101. We add 
additional contributions to the reported uncertainty, in 
quadrature, to allow for the absolute flux calibration 
and the colour corre ction uncertainty as described in 
iSevmour et al.l (|2007l ). We supplement this photometry 
with new observation s from Hersch el using the PACS 
(jPoglitsch et al.l and SPIRE (IGriffin et al.l 

instruments, as well as a LABOCA/870 /im flux density 
(Dannerbauer et al., in prep., ESO program ID 084. A- 
1016(A) and MPG program ID 083.F-0022 - both PI 
J. Kurk). The Herschel observations were carried out as 
part of the project scientist guaranteed time (PI Altieri). 

The SPIRE observation (OBSID 1342210877) was per- 
formed using a four repetition cross-scan Large Map 
mode covering 10 x lOarcmin, achieving 1 a instrumen- 
tal noise of (3.5, 4.8, 5.4) mJy at (250, 250, 500 /im) - 
slightly below the S PIRE extragalactic confusion noise 
([Nguyen et al.ll2010D . To obtain SPIRE photometry for 
the Spiderweb Galaxy we used a time-line fitteiEJ which 
provides a better estimate of the flux density, because no 
binning into sky pixels is involved in the process. The 
final flux density uncertainty, reported in Table [IJ in- 
cludes the fit error, the confusion noise and a conserva- 
tive 7 % flux calibration uncertainty (see the SPIRE Ob- 
servers' Manual). The photometry from standard source 
fitting on maps is consistent with the time-line fitter re- 
sults, but we prefer to use the time-line fitter results as 
these are generally more accurate, especially for sources 
at high signal-to-noise. The PACS observations (OBSID 
1342222456 and 1342222457) consisted of two crossed 
scan maps covering the same area and the data were re- 
duced using HIPE v6.0 (|Ottll2010D . The data cubes were 
processed with a standard pipeline, including a sliding 
high-pass filtering on the detector time-lines to remove 
detector drifts and 1/f noise with an iterative masking 
of the sources. We extracted flux densities via stan- 
dard aperture photometry techniques and the PACS un- 
certainties include contributions from instrumental noise 
and absolute flux calibration. For the LABOCA/870 /mi 
flux density uncertainty we add a 15 % uncertainty, in 
quadrature, to the instrumental uncertainty to account 
for absolute flux calibration, atmospheric effects. 

The Herschel beam varies from 6 — 36 arc-seconds, but 
the radio galaxy appears point like in all-bands from 8 /im 
to 250 /im (see Fig. p. The 350, 500 and 870 /on im- 
ages all show a slight extension in the east-west direc- 
tion consistent with emission from the clump of galax- 
ies approximately 30 arc-seconds to the west of the radio 
galaxy seen in the 8 /im through 250 /im bands. This ex- 
te nsion is consist e nt wit h the SCUBA/850 /im detection 
bv iStevens et al.l ([2003). This clump comprises of four 
sources seen individually at 8 /im and 24 /im, the bright- 
est two of which have been spectroscopically confirmed 
to lie at t he same redshift a nd therefore be proto-cluster 
members (jKurk et al.l 12004?) ; hence, we can assume the 
long-wavelength flux density from these sources to have 
similar colours due to the Rayleigh- Jeans tail. We note 



the peak flux density at the longer wavelengths remains 
centered on the radio galaxy. We assume, for the pur- 
poses of SED fitting, that all this flux density comes from 
one source, but then treat the SFR measured from the 
starburst component as an upper limit to that from just 
the radio galaxy. 

All the data used in the fitting of this sample are pre- 
sented in Table [TJ 

3. SED FITTING 

The goal of our SED fitting is to obtain the AGN and 
starburst IR luminosities (8 — 1000 /im) and then derive 
physical quantities from these results. We fit the ob- 
served photometry across 5 — 300 /im rest frame with 
a two-component model consisting of an AGN and a 
starburst template. Extrapolation of the radio core flux 
density using the observed flux densities at wavelengths 
>6cm and the core spectral index from Caril li et al.l 
(|1997|) . a = —1.3 (S u oc v a ), suggests that synchrotron 
emission at these wavelengths is many orders of mag- 
nitude fai nter. We use the Deco mpiiPl SED model fit- 
ting code (jMullanev et al.l 120111 ) to fit the two compo- 
nents with some minor adaptations to use additional 
band-pass filters. The AGN component is an empiri- 
cal model based on observations of moderate-luminosity 



2-inkBV 



1 4 ' 2 - lO^ergs" 1 ) low-redshift AGN 
in iMullanev et al.l (|2011[ ) selected at high X-ray energies 
from the local, z < 0.1, Swift-B AT sample (jTueller et al.1 
2008). This model is a continuous function comprised of 
two power-laws and a modified black-body, expressed as: 



F„ oc 




6 /im < A < 19 /im 
19 /im < A < 40 /mi 
A > 40 /im 



(1) 



In comparison with other AGN templates commonly 
used, it is broadly similar in the mid-IR regime; however 
it has lower far-IR emission, when normalized at 19 /im, 
compared to NGC 1068 and Mrk 231, likely attributable 
to unaccounted starburst contri butions at far-IR wave - 
lengths in these two galaxies (|Le Floc'h et al.l [2001). 
This model is also broadly similar t o e mpirical mod- 
els pre sented by iRichards et al.1 (|2006[ ) and iNetzer et al.1 
(2007). The advantage of the formulation above is that 
the various parameters (the two power-law slopes and the 
black-body peak wavelength) can be adjusted to closely 
match these other models which we test later on. 

The starburst models consist of five SED s developed 
by Mullaney et al. from lBrandl et ail (|2006[ ) designed to 
represent a typical range of SED types (i.e., peaking at 
a range of wavelengths). Beyond 100 /im these templates 
are extrapolated using a grey body with a frequency de- 
pendent emissivity approximated as ~ v$ with ft = 1.5. 
In addition to the normalisation and choice of starburst 
template, Decompir can apply extinction to each com- 
ponent separately. Asthe AGN-dominatcd mid-IR spec- 
trum (Og le et al.l l2012f ) shows a little silicate emission, 
not absorption, we only allow the starburst component 
to be affected by absorption^ Therefore the chi-squared, 

24 https : / /sites . google . com/site/decompir/ 

25 We note that while the X-ray emission fro m AGN shows mod - 
erate extinction, N H =(2.6 ± .0.6) X 10 22 cm~ 2 llCarilli et al.i r2002), 

herschel. esac. esa. int/twiki/bin/view/Public/SpireCallbrationWeb the mid-IR SED does not IQgle et al J 120121 ) . This is consistent 
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\ 2 , minimisation performed by Decompir has four free 
components: normalisation of each template, choice of 
starburst template and starburst extinction. 

The best fit combination of models returns a reduced 
chi-squared of Xr = 1-4 with the starburst template, 
'SB4', the template peaking at the shortest wavelength. 
We note that while the best fit model has some extinction 
for the starburst component, it is poorly constrained and 
we get an equally good fit with no extinction and similar 
luminosities within 2%. The fit is shown in Fig. [2] and 
the luminosities are reported in Table [2] The model fits 
the data well, although at the longest wavelengths the 
model slightly struggles to fit the sharp decrease in flux 
density from 500 to 870 fim. The SED at these wave- 
lengths is dominated by the Rayleigh- Jeans tail of the 
starburst component and any blackbody would struggle 
to fit such a flux density ratio, although the slope does 
depend on the emissivity index. The 'SB4' starburst tem- 
plate is the hottest of the five used, with a dust temper- 
ature of Tdust ~ 64 K. The next two cooler templates, 
Tdust ~ 46 - 53 K, give a xl ~ 5 - 8 and fit the 500 fim 
flux density better, but exceed the 870 fim. flux density 
by > 4cr. 

If we exclude the 500 jtim data point we marginally 
improve the best fit to Xr — 1-1, again for the 'SB4' 
template and with the output parameters changing very 
little (note the number of degrees of freedom drops by 
one). One potential method to better fit the observed 
high 500 to 870 fim flux ratio is to use starburst SEDs 
with a steeper emissivity index. By modifying the tem- 
plates so that (3 = 2 we can obtain a steeper slope on 
the Rayleigh- Jeans side of the SED peak. Repeating the 
SED fitting with these modified starburst templates we 
obtain a best fit of Xr = 1>1 from the modified 'SB5' 
starburst template. The overall fit does not significantly 
change with a slightly smaller 500 /im excess, by ~ 20 %, 
as well as an AGN luminosity increase of ~ 10 % and a 
starburst luminosity decrease of ~ 20 %. Hence, there is 
no significant improvement with this change to the star- 
burst templates. 

Visual examination of the images at 350, 500 and 
870 jtim, see Fig. Q] suggests the measured flux densities 
at these wavelengths may have contributions from other 
proto-cluster members as discussed in iJ5J As the galaxies 
responsible for this emission lie at the same redshift, and 
are almost certainly starbursts to be emitting at these 
wavelengths, they likely have similar flux density ratios. 
The 250 /im flux density is not affected by these nearby 
sources and the peak flux density still originates from 
the central radio galaxy. The ratio of the radio galaxy to 
companion peak flux densities at 250 /zm is about 2:1; 
hence, we estimate the contribution of the companion to 
the measured flux density at longer wavelengths to be 
~ 1 at most. We note that the 350 /jm flux density need 
only be decreased by 25 % to produce an unphysical flux 
density ratio for a cold dust SED. However, if we decrease 
the three longest wavelength flux densities by this much 
the best fit gives a starburst luminosity only 10 % lower. 
Hence, although the estimated starburst IR luminosity is 
formally an upper limit for the radio galaxy, we postulate 

with the bright H a seen by Ncsvadba ct al. ( 2006|) and can be at- 
tributed to differences in dust grains seen toward the circumnuclear 
regions of AGN and the diffuse ISM of our galaxy. 



that the true starburst IR luminosity of the Spiderweb 
Galaxy is, at most, only 10% lower. 

We do not have enough photometry at 5 fim < A rcs t < 
50 /im to allow all the AGN parameters (the two power- 
law slopes and the black-body peak) , other than the nor- 
malisation, to vary simultaneously. Allowing just one 
parameter to vary results in changes of luminosities of 
< 5 %. W e can change the param ete rs in equation one to 
mimic the lRichards et al.1 (|2006[ ) and lNetzer et all (2007) 
AGN models. In both cases, the AGN luminosity halves 
and the starburst luminosity doubles, but the fit is worse 
with xl ~ 8-5. We note that this fit is improved if we use 
the Richards or Netzer AGN model with the /3 = 2 star- 
burst SED, with the same qualitative results, but with 
a Xr = 3.3. We also trie d using the s uite o f empirical 
starburst templates from iRieke et al.1 (|2009| ) instead of 
the Brandl templates. The best fit templates were also 
those peaking at the shortest wavelengths, i.e., the ul- 
tra luminous IR galaxy (ULIRG) templates, but with a 
slightly higher chi-squared: Xr ~ 3. To obtain a better 
understanding of the uncertainties we also randomly var- 
ied the flux densities with a Gaussian distribution based 
on their quoted uncertainties one hundred times and re- 
peated the fitting procedures. We found the variance of 
the resulting luminosities to be < 10 %, i.e., less than the 
quoted uncertainties in Table [5J 

We finish this section by noting that while there are 
many ways to tweak both the AGN and starburst tem- 
plates and change the parameters of the fit all solu- 
tions have an AGN component dominating at the short- 
est wavelengths, A ros t < 30 /im, and a starburst com- 
ponent dominating at long wavelengths A rGS t > 60 /im. 
Hence, all good fits have a powerful, Ltr(SB/AGN)> 
5 x 10 12 Lq, starburst and AGN component. Until we get 
high spatial and/or spectral resolution data, e.g. from 
ALMA or JWST, it will be difficult to differentiate be- 
tween many moderately good fits. 

4. DISCUSSION 

We have measured very high AGN and starburst lumi- 
nosities which imply a high central black hole accretion 
rate and SFR, respectively. The total IR luminosity of 
this system causes it to be classed as a hyper-luminous 
IR galaxy, Lir > 10 13 Lq . We discuss here what these 
results mean and how they relate to the radio emission 
and overall evolutionary state of the system. We also 
discuss the potential origin of the 500 /jm excess. 

4.1. AGN Power 

The total IR AGN luminosity marks the Spiderweb 
Galaxy as one of the most powerful AGN known; there- 
fore the AGN is likely having a strong effect on its 
host galaxy and immediate environment (see ^4.6|) . As- 
suming the full u nobscured AGN SED is similar to the 
lElvis et all (|1994f ) templates, we get L BO l(AGN) - 6 x 
Lir(AGN) (we note other unobscured AGN templates 
produce similar numbers) and therefore an AGN bolo- 
metric luminosity of ~ 7x 10 13 Lq. This bolometric AGN 
lumin osity estimate is high er than previous estimates, 
e.g., iNesvadba et al.l (120061) who used the observed X- 
ray lu minosity (iCarilli et al.l 120021 ) and corrected it us- 
ing the lElvis et al.l (|1994[ ) template. L B ol(AGN) - 10 x 
£ 2 -iokev(AGN). 
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This discrepancy is unlikely due to extinction as seen 
towards the nucleus in X-r ays given the low absorbing 
column (Caril lTet al.l I2002T) . As noted earlier, neither 
the mid-IR spectrum (fOgle et al.l 12012ft nor the SED 
(Drouart et al. 2012) indicate significant extinction ei- 
ther. The discrepancy between the hard X-ray and IR- 
derived bolometric AGN luminosities is more likely due 
to the dispersion between individual sources and the 
lElvis et aLl (|1994l) templates. This, in turn, could be 
explained by a different dust-to-gas ratio in this source, 
a clumpy torus or an unusual AGN geometry, but further 
investigation is beyond the scope of this paper. 

This galaxy exhibits both a large scale central ve- 
locity distribution in dicative of a powerful outflow 
(|Nesvadba et al.ll2006l ) a nd a high-lumino sity mid-IR H 2 
molecular emission line (|Qgle et al.ll2012h . The derived 
bolometric radiative power would theref ore be enough to 
powe r both the outflow (~8x 10 12 L Q - INesvadba et al.l 
120061) and the prob ably shock-heated gas (~4x 1O 1O L 
- IQgle et al.l I2012T) . even assuming a 10 % coupling ef- 
ficiency. However, the kinetic power of the radio jet, 
estima ted to be 0.6 — 8.0 x 10 13 L Q by INesvadba et al.l 
(|2011|) . is also powerful enough to potentially drive the 
outflow. 

Nesvadba et al. (2006) also present an estimate of 
the black hole mass from the measurement of the central 
H a luminosity and FWHM. The principle uncertainty 
in this measurement is how th e orientation of the g alaxy 
affects the application of the iGreene fc Hoi (|2005f) scal- 
ing relationship, which is based on a sample of lower 
luminosity AGN likely to be more face-on than the Spi- 
derweb galaxy, i.e., at lower inclinations. Indeed, G. 
Drouart et al. (2012, in preparation) find a correlation 
between radio core prominence and viewing angle from 
a hot torus model fitted to the mid-IR SED of HzRGs 
which suggests that the radio axis is aligned to torus 
axis. These authors find that the Spiderweb Galaxy is 
inclined at angle of ~ 68deg to the line of sight, but 
has modest extinction, an Ay ~ 1.5. While generally 
there may be misalignments of up to 10 — 20deg be- 
tween the radio and torus axes, it is likely that the Spi- 
derwe b Galaxy is more edge-on than th e IGreene fc Hoi 
(|2005D sample. IMcLure fc Dunlopl (|2002f ) show that the 
measured FWHM from a thin accretion disk is propor- 
tional to the sine of the inclination angle squared. There- 
fore, if the Spiderweb Galaxy has a higher inclination 
angle, 0rg, than the mean value for the AGN sample 
used in IGreene fc Hoi ((2005), 6*gho5j then we would be 

overestimating the black hole mass by sl s " n i^ ™^ which 

INesvadba et al.l (|2011[) assume to be ~ 2. We there- 
fore take the lower, inclination-corrected, value of the 
black hole mass from lNesvadba et al.l IpOllh . ~ 10 10 M . 
The higher bolometric luminosity derived above implies 
a higher black hole mass derived from the H a lumi- 
nosity. However, we also note the black hole mass re- 
ported here would put this galaxy on par with some of 
the m ost massive black hole s known in the local Uni- 
verse (jMcConnell et al.|[20TTl ). and are very rare for the 
z > 1 Universe. Hence, we treat this black hole mass as 
uncertain by a factor 2 — 3. 

We then relate the bolometric luminosity to the max- 
imally possible accretion luminosity where the contin- 
uum radiation force outwards balances the gravitational 



force inwards in hydrostatic equilibrium, i.e., the Ed- 
dington luminosity, and we find an Eddington accre- 
tion rate of ~ 20 %. This is a factor of a few higher 
than the estimate of INesvadba et all (j 2 1 lh due to the 
higher bolometric luminosity and lower black hole mass 
used in this work. This high rate of accretion is com- 
mon in less bolomet rically luminous AGN, (~ 25 % - 
iKollmeier etall l2006) . but is rare for less radio-luminous 
AGN (Lsoomhz ^ lO^WHz" 1 - IBest fc Hecln7a^l27)12h 
at least in the local, z < 0.1, Universe. However, it 
is p robably common for HzRGs as b ased on their X- 
rav (lOverzier et al.N2005l ) and mid-IR (|Qgle et al.l 120061: 
ISeymour et al.ll2007l ) luminosities. These results suggest 
the accretion disk is in an efficient state, although the 
large outflow hints that it may soon run out of fuel. 

4.2. Star- Formation Rate 

The starburst IR lu minosity implies a SFR of 1390 ± 
150M Q yr _1 using the iKennicuttl (|1998f ) scaling relation 
(and therefore a Salpeter IMF). We noted in f}2] that 
the longer wavelength bands may have some contribu- 
tion from nearby proto-cluster members, but estimate 
that this contribution is < 10% in P This SFR is 
much higher than estima tes from optical/near- IR data 
(|Hatch et alj 120081 120091 ), likely due to obscuration by 
dust, but is broadly consistent with estimates from the 
PAH luminosities measured from the mid-IR SED (Ogle 
et al., 2012; J. Rawlings, 2012, in preparation). We note 
that the IR luminosity is more sensitive to the star for- 
mation averaged over a slightly longer time frame than 
optical/UV observations. We can estimate the mechani- 
cal energy o f this starburst using the canoni cal conversion 
factor from lDalla Vecchia fc Scha ve (2008;) and obtain a 
value 1 — 2 x 10 11 L© which is ins ufficient to explain th e 
kinematics and outflow observed (Ncsvadba~e£jd][2006) . 
However, any material driven out may fall back onto the 
central galaxy or the satellite galaxies over time (e.g., 
lHatch et al.ll2009l and see p~6]) . 

Estimating the stellar mass of this so urce is diffi- 
cult d ue to the strong AGN continuum. lHatch et al.l 
( 2009) provide an estimate of 1.1 ± 0.2 x 10 12 M Q 
using SED fitting of high resolution HST rest-frame 
UV/optical data. This result is consistent with the 
upper limit found from rest-frame near-IR o bservations 
(ISeymour et al.l|2007t IDe Breuck et aLll2010l) . We adopt 
the lHatch et al.l (|2009f ) value as an upper limit since it 
may still be contaminated by the AGN as it lies in the 
AGN-do minated region of mid-IR colour-colour space as 
shown in IDe Breuck et all 1)20101 ). We can now derive 
the specific SFR, sSFR=SFR/M 5te iiar, and find a value 
of > 1.3 Gyr -1 . This is a high sSF R, but common for 
star forming galax ies at this redshift (Noeske et al.ll2007l : 
lElbaz et al.l 1201 J ) . However, it is unusual for a galaxy 
this massive. 

The inverse of the sSFR gives the characteristic time, 
^char, in which the galaxy would double its stellar mass 
if it maintains the current SFR. While this time-frame is 
short, < 0.8 Gyr, for the Spiderweb, it is likely that the 
star formation will turn off on even shorter timescales. 
Measurements of CO provide the best estimates of the 
mass of molecular hydrogen, i.e., the main component of 
the fuel for star for mation and black hole growth (e.g., 
lEmonts et al .1 120111 ) . While such measurements do not 
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currently exist for the Spiderweb Galaxy we infer a value 
of Mh 2 ~ 0.5 - 5 x 10 10 M Q based the far-IR lum i nosity 
and the conversion factors used in lEmonts et al.l ([20111 ) 
assuming the Spiderweb is similar to the HzRGs in that 
study. If we assume 100 % star formation efficiency, then 
this SFR can only be maintained for another 4 — 40 Myr, 
neglecting infall from satellite galaxies which may hap- 
pen within a few hundred Myr (see below). If the mass 
of material esti mated to being remov ed by the outflow, 
> 4OOM yr _1 ([Nesvadba et al.ll2006D . is taken into ac- 
count, this time frame decreases by ~ 30%, more if the 
powerful jets also entrain and remove significant amounts 
of gas from the galaxy, or if the molecular hydrogen mass 
is lower than assumed. If we assume a much lower star 
formation efficiency, e.g. 10%, then this time-scale de- 
creases by an order of magnitude, but leaves fuel for 
further bursts of star-formation triggered by its proto- 
cluster environment. 

4.3. Age of Radio Jet 

Previously. INesvadba et al] (|2006f ) estimate d the age of 
the ra dio jet using the method put forth by IWan et all 
(2000) and obtained a value of 4 - 10 Myr. We can 
re-estimate the age us i ng mo dels of radio lobe evolu- 
tion from IMocz et all (|2011l ). These authors repro- 
duce the emission at radio wavelengths from synchrotron 
processes and the X-rays emitted by the lobes due to 
inverse Compton (IC) scattering of the Cosmic Mi- 
crowave Background. These models lead to tracks of lobe 
radio/151 MHz and X-ray/ 1 keV luminosity against time 
and si ze of radio lobes at different redshifts. ICarilli et al.l 
(2002) estimate that 10 — 25% of the X-ray emission 
they measure is spatially extended on scales of 10 — 
20 arcsec, similar to the extent of the radio jet, 15.8 arcsec 
(|Carilli et al] 119971) . Hence, if we assume all this X- 
ray luminosity is from IC processes we get a lobe X- 
ray luminosity of L2-iokeV ~ 1 — 2.5 x 10 n L Q . As 
IMocz et all (f20lTh use a 1 keV X-ray luminosity wc take 
LikcV ~ 1-8 — 4.5 x 10 11 Lq (assuming a photon power- 
law index of T = 1.4), and from the 365 and 74 MHz ra dio 
flux densities (|Douglas et alj|1996tlCohen et al]|2007D we 
find a 151 MHz luminosity of L 151 = 5.6 x 10 29 WHz" 1 - 
2.2 x 10 11 Lq. Using the redshift and inclination estimate 
from G. Drouart et al. (2012, in preparation) we estimate 
the total lobe size to be ~ 140 kpc , which implies an age 
of - 5 Myr in Fig. 1 of lMocz et all (|201lD using the more 
powerful, z — 2 AGN track. Using the radio and X-ray 
luminosities, where L151 > LikcVi m Fig- 3 of Mocz et 
al., we obtain an estimate of < 25 Myr for the age of the 
lobe. Hence, we take the age of the lobes to be in the 
range 5 — 25 Myr. 

4.4. Possible [CII]158fj,m Line 

In SJ3] we reported an excess of emission in the 500 /im 
band which we measure to be 10.0 mJy at > 1 a signif- 
icance (~ 2cr based on just the instrumental noise). A 
potential explanation for this slight excess could be a 
significant contribution from the redshifted [CII]158/im 
line (A bs = 499 /im). We note that no oth er lines, e.g. 
CO, f all within the 500 /im band width. ISmail et al] 
<|201 If ) examine the effect of strong far-IR emission lines 
on broad band photometry and show that such emis- 
sion lines can cause measurable excesses in the broad 



band photometry. The excess measured here, 40% above 
the 25.3 mJy of the model, corresp onds to a value o f 
L\cn]/Lm ~ 2.7% (e.g. Fig. 4 of ISmail et al.ll20lll) . 
While this value is higher, by a factor of ~ 2, than any 
known value, it potentially indicates evidence for strong 
[CH]158/im distributed across the proto-cluster environ- 
ment, due to the extended 350 — 870 /im emission in 
the direction of confirmed proto-cluster members. The 
idea of a high, extended [CII] luminosity is consistent 
with the ide a of a gaseous en vironment of neutral gas, 
ioniz ed gas (jKurk et al.l l2003h . shock ed Ho (lOgle et al.l 
l2012f ) and intra-galaxy star formation ([Hatch et a,l.ll2009f ) 
within this proto-cluster environment. We caution that 
the statistical significance of this result is poor but the 
Spiderweb Galaxy remains a promising target for far- 
IR/millimetre spectroscopy. 

4.5. Evolutionary State of the Spiderweb Galaxy 

We have determined that the Spiderweb Galaxy is in 
a rare, short-lived, phase of extreme, concurrent black 
hole and host galaxy growth. If we compare the relative 
luminosities of the AGN and sta rburst component to the 
'continuous growth' relation in iLutz et all ((20l3), i.e., 
the ratio required to maintain a constant black hole to 
stellar mass ratio, we see that the black hole is growing 
at a relatively faster rate by almost 1 dex. This result 
is consistent with the observed, albeit rather uncertain, 
ratio of stellar to black hole mass, < 110, i.e. the black 
hole is relatively more massive (and has been growing 
rapidly) than it would be for the Spider web Galaxy to lie 
on th e local i M— a' mass ratio of ~ 200 (jMagorria n et al] 
1998). The stellar an d black hole masses do lie on the 
relationship found by lHaring fc"R ix (2003) and we note 
that both observed relations have a large scatter, 0.3 dex. 
If this source is to end up on the local 'M — a' relation, 
further build-up of stellar mass is required. However, we 
have estimated that the star formation rate will likely 
rapidly decrease on time-scales of ~ 40 Myr assuming no 
f urther sourc e of fuel . 

IPage et all (|2012f ) have recently examined the 
Herschel/250 /im detection rate of X-ray luminous AGN 
in the Chandra Deep Field North in the 1 < z < 3 red- 
shift range as a function of X-ray luminosity. Those au- 
thors found that the fraction of X-ray AGN detected with 
"5250/im > 18mJy decreases to zero above L2_8koV = 
10 44 ergs _1 (= 3x10 10 Lq) and interpreted this result as 
suppression of star formation by powerful outflows driven 
by the AGN. The Spiderweb is a clear counter-example of 
this result with an X-ray luminosity well above this limit 
and a brigh t 250 u rn detection, S250 = 42 mJy. However, 
IPage et all ([20121 ) surveyed a small area, covering just 
0.12 deg 2 , and hence only probed a small volume where 
galaxies as massive as the Spider web are rare over this 
volume (e.g., IFontana et al"1l2006D . We can estimate an 
upper limit to the duty cycle of a simultaneous AGN 
and starburst activity above these limits by dividing the 
look-back time probed across 1 < z < 3, 3.64 Gyr, by 
the number of X-ray AGN in this study, 21, and obtain 
a value of 170 Myr. If the Spiderweb Galaxy is in a tran- 
sition from merger-induced starburst phase to an AGN 
phase, and this transition is relatively short, < 170 Myr , 
as the analysis here suggests, then the IPage et all ([2012D 
results are consistent with our observations. We point 
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out that such comparison is only valid if there are sim- 
ilar evolutionary processes going on in the less massive 
iPage et al.l (j20i2fl sample. 

4.6. Proto-Cluster Environment 

The environment of the Spiderweb Galaxy has been ex- 
tensively studied and a full review of work in this area is 
beyond the scope of this paper. We have noted previously 
that this galaxy hosts one of the most powerful AGN 
known: both from its bolometric luminosity (indicative 
of radiative power) and its radio lumin osity (indicative 
of kinetic power). iTanaka et al.l ([2010D compared mas- 
sive galaxies in the environment of the Spiderweb Galaxy 
with identically selected field galaxies from the GOODS- 
S field. They determined that the galaxies in the proto- 
cluster environment had formed earlier and had lower 
SFRs, particularly close to the HzRG and along the east- 
west distribution of proto-cluster members. The radio 
emission has the same projected east-west axes; hence, 
it is possible that the lower SFRs observed could be due 
to the influence of the AGN jet. Alternatively, they may 
simply be more heavily obscured. 

As the Spiderweb Galaxy is situated in a well studied 
proto-cluster we can estimate its future growth. Sev- 
eral authors have identified proto-cluster members and 
estimated that they would merge with, or be tidally 
stripped by, the ra dio galaxy before the present epoch. 
lHatch eta l. (2009) identified 18 proto-cluster members 
through HST imaging within 150 kpc of the central ra- 
dio galaxy, the combined mass of which they estimated to 
be 1.9±o | x 10 11 M e . iDohertv etahl (|2010f ) spectroscop- 
ically confirmed two red galaxies within ~ 500 kpcPl. in 
addition to those mentioned above, with stellar masses 
of the order 3 — 5 x 10 11 M© each. Hence, the total 
stellar mass increase of the Spiderweb galaxy over the 
upcoming ~ 10 Gyr would be ~ 10 12 M©, i.e., at least a 
doubling in its mass, which would put it closer to the lo- 
cal l M — a' relation. iKuiper et al.l ([20111 ) interpret their 
kinematic study of the environment, which showed a dou- 
ble peaked velocity distribution, as being due the merger 
of two massive clusters. If each cluster contained a dom- 
inant central galaxy then the high AGN and starburst 
luminosities could well be due to a gas rich merger of 
two large galaxies. 

5. CONCLUSIONS 

We have measured the total IR rest-frame SED of the 
Spiderweb Galaxy, for the first time, and fitted an AGN 
and starburst model components. We find that the bolo- 
metric luminosity of the AGN component of this source 
implies that the AGN has an accretion rate high enough 
to power both the jets and the observed outflow. The 
Spiderweb Galaxy has a high SFR; however, we estimate 
this SFR cannot be maintained for long and cannot be 
responsible for the outflow. We estimate the radio jets 
are 5 — 25 Myr old. At the current SFR, the host galaxy 
has only had enough time to accumulate ~ 6 % of its 
observed stellar mass while the jets were on. Hence, the 
jets were likely triggered after most of the current burst 
of star formation commenced. The central black hole 
is possibly growing at a relatively faster rate than the 

26 Not the ~ 300 kpc previously reported. 



host galaxy (assuming this system will end up on the 
local M — a relation). Studies of the proto-cluster en- 
vironment suggest that mergers of satellites, and hence 
increased star formation, will occur on timescales of hun- 
dreds of Myr and ultimately the stellar mass will double 
by the present epoch. We also see very tentative evidence 
for strong [CII]158/im emission. 
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Table 1 

Photometric data 



wavelength flux density reference 

fim mJy 

~~24 3.89±0.25 (I) 

100 28.5±3.1 (2) 

160 43.0±2.9 (2) 

250 42.0±7.2 (2) 

350 38.9±7.9 (2) 

500 35.3±8.3 (2) 

870 6.7±1.3 (3) 



References. — (1) Dc Brcuck ct al. 
{HUD; (2) this work ; (3) Dannerbauer 
et al. (2012, in preparation). 



Table 2 

Derived physical properties 



property 


value 


reference 


L IR (Total) 


1.97 ±0.28 x 10 13 L 


(1) 


L IR (AGN) 


1.17 ±0.27 x 1O 13 L 


(1) 


L IR (SB) 


0.79 ± 0.09 x 10 13 Lq 


(1) 


Af gal (stellar) 


< 1.1 x lO i2 M 


(2) 




~ 1 X 1O 1O M 


(3) 


SFR (optical) 


57±8M yr- 1 


(2) 


Inclination, 8 


< 68 deg 


(/D 


Angular Size 


15.8 arcsec 


(o) 


L 2 -iokcV (unob.) 


< 1.0 x 10 12 L Q 


(6) 


L BO l(AGN) 


~ 7.0 X 10 13 Lq 


(1) 


^Edd 


~ 0.2 


(1) 


SFR 


1390 ±150 Mgyr" 1 


(1) 


sSFR 


> 1.3 Gyr" 1 


(1) 


^char 


< 0.8 Gyr 


(1) 


D/sin(6) 


~ 140 kpc 


(1) 


Ll51MHz 


5.6 x 10 29 WHz" 1 


(1) 


^lkcv (lobe) 


1.8-4.5 x lO n L 


(1) 


tjet 


5-25 Myr 


(1) 



References. — (1) this work; f2l lHatch eFaTI l[2009h 
— note the SFR is int egrated over a 32 kpc radius; (3) 
INesvadba eTall ll20Tll); (4) Drouart et a l. (2012, in 
prepa ration); (5) Carilli ct al. (1997); (61 ICarilli eTall 
J20021 ). 

Note. — Results of SED fitting are listed first, pa- 
rameters from the literature are listed second, and de- 
rived measurements in this work are listed last. Total 
IR luminosities are integrated over 8 — 1000 fira. SFR 
is based on the Kennicutt (1998) conversion which uses 
a Salpeter IMF. 
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Figure 1. Postage-stamp IR images, 70 X 70 arc-seconds, of the Spiderweb Galaxy (RA 11:40:48.60, Dec. -26:29:8.50; J2000) with North 
up and East to the left. From top left to bottom right the images are: Spitzer/ 8/ 16/24 /im, PACS/100 and 160 /mi, SPIRE 250/350/500 /mi 
and LABOCA/870/tm. The 8 and 16 /mi band are shown for illustration only and are not used in the model fitting. Note the Airy ring is 
quite prominent in the Spitzer data. 
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Figure 2. SED fit to the rest-frame 5 — 300 /im rest-frame photometry (black dots) of the Spiderweb Galaxy at z = 2.156 (MRC1138-262). 
The total model fit, and AGN and starburst components are indicated in the figure. The vertical grey dotted line represents the position 
of the [CII]158/jm line. 



